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Abstract 



The methods of neutrino spectra calculations with GEANT library and the peculiarities 
of such calculations for different neutrino beams are summarized. The formula expressing 
K*^ yields from a neutrino target through the if + and K" yields from the same target is 
proposed. 
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1 Introduction 



In 1988 - 1989 the new generation of neutrino experiments started at IHEP 
(Protvino). It was a result of the physical start of the new IHEP-JINR neutrino de- 
tector, increased intensity of the IHEP 70 - GeV proton accelerator and somewhat 
changed situation in the neutrino physics, which required better accuracy of mea- 
surements and inspired interest to experiments with electron neutrinos. 

Most of these experiments were carried out with a nonstandard mode of work of 
the neutrino facility. It concerns the beam-dump experiment, the experiment with 
focusing system switched off and the new experiment with short decay base. These 
experiments required a new program for neutrino spectra calculations which could 
take into account additional effects that become important in these modes of the 
neutrino beam formation, for example hadronic interactions in the elements of the 
neutrino channel other that the neutrino target. 

For this purpose a program based on the GEANT library [1] was chosen. It has 
the following advantages: 

• developped, multiply checked and supported package for a simulation of passage 
of particles through matter, including hadron showers; 

• the way of describing a geometry is sufficiently developped and familiar to most 
of physicists. 

The disadvantage of the full GEANT - like simulations of processes in the 
neutrino channel is that it is rather time consuming. This required usage of tricks 
that would increase the number of useful decays, but in spite of them, only after 
appearance of sufficiently powerful computers an accurate calculation of neutrino 
spectra by such a program became possible. 



2 The program for the neutrino spectra calculation 

In the GEANT library the procedures for random decays of particles in flight are 
provided. If GEANT accomplishes a decay, the particle is stopped (no more traced). 
Thus, a particle can decay only once, which required the simulation of many parti- 
cles and a lot of CPU time. For this reason routines different from the ones of the 
GEANT library were used for the simulation of decays. These routines are called 
several times on every step of a particle. The GEANT table of branching ratios was 
changed (some branching ratios are zeroed), so that GEANT does not accomplish 
decays to neutrinos of the type for which the spectra are being calculated. To take 
into account correctly the zeroed branching ratios the GEANT routine GDECAY 
is slightly changed. Namely, if no decay mode is chosen (it is now possible since 
the sum of branching ratios is not equal to unity), a particle is not stopped. Thus, 
a particle makes more steps, which also improves the statistical accuracy of the 
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calculations. The effect of decreasing number of particles along the beam due to a 
finite lifetime is taken into account through weights ascribed to particles. 

Unfortunately, although all these tricks certainly decrease the overall statistical 
errors of calculations, they do not decrease by the same factor all components of 
these errors. For example, the fluctuations due to hadron showers are not decreased. 
This makes very difficult a direct estimation of statistical uncertainties of spectra. 
The only obvious way to estimate them is to compare results of several runs of the 
program. This consideration also suggests that there is no sense in very small steps 
or many decays in one step. 

Some branching ratios, such as the branching ratios to hadrons, are not changed 
in order to simulate more exactly the beam composition (for example if a kaon de- 
cays to pions, these pions are also traced later). For the same purpose and in order 
to take into account decays of muons during the calculation of electron neutrino 
spectra the program has 2 modes of work: 

(1) The calculation of muon neutrino and antineutrino spectra. The branching ra- 
tios of decays to electron neutrinos are not changed since in some of them 

charged pions are produced. 

(2) The calculation of electron neutrino and antineutrino spectra. The branching 
ratios of decays to muons are not zeroed. 

In each of these modes it is possible to simulate directly hadron showers in the 
neutrino target induced by initial protons or (not for the beam dump) to use param- 
eterized yields of particles from the standard neutrino targets. 

As most of programs using GEANT library the program consists of the follow- 
ing parts: 

(1) The initialization. Here beside usual GEANT calls the GEANT table of branch- 
ing ratios is changed. 

(2) The geometry definition. It is a usual GEANT geometry definition [1] except 
that some volumes should have standard numbers (the target, lenses and some 
volumes for the calculation of control distributions). Besides, the constants 
for the spectra calculations, such as a distance to a possible detector and its 
maximal sizes, are set. 

(3) The routine for the action before each event. Here either a proton is put into 
the GEANT initial bank (direct calculation) or several hadrons with random 
momentum and angle (parameterization). 

(4) The routine for the action after each step. It calls the two main routines of the 
program: 

a) The routine that ascribes a weight to a new particle or calculates a weight 
after each step. A weight 1 is ascribed to primary protons. In the mode of work 
with a parameterization a weight to primary hadrons at their first appearance 

is ascribed according to this parameterization. For other particles, at their first 
appearance, their creation vertex is searched for in the array of vertices (this is 
not in the standard GEANT banks and common blocks, see below). From this 
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array the initial weight, the type and momentum of a parent particle, the spin 
direction (for a muon) are taken. If there is a correction factor for the process 
of creation of this particle by its parent particle the weight is corrected. After 
each step the weight is decreased according to the life time of particle divided 
by the sum of all zeroed branching ratios. 

b) The calculation of the neutrino spectra. For two body decays the limits on 
the angles in the center of mass system are calculated, so that neutrinos from 
a decay can be inside the maximal detector volume. For three body decays, 
if a neutrino energy is above threshold and its direction does not hit the de- 
tector maximal volume, its azimuthal angle is resimulated withing the limits 
that contain the maximal detector volume. The probability of a decay in case 
of additional limits is multiplied by the corresponding phase space factor. The 
matrix element of decays is taken into account. For decays of muons it de- 
pends on the spin direction. 

The probabilities of decays are put into double precision arrays that corre- 
spond and will be rewritten in the end of run into one- and two-dimensional 
histograms (energy - radius, energy - angle). 
(5) The routine for the action after each act of creation of new particles. As usual, 
particles below cut and those that do not contribute to the neutrino spectra, are 
ignored. Here the array of vertices mentioned above is filled 



3 The spectra of hadrons. 



The first calculations of WBB spectra were made in the mode with the full simu- 
lation of the hadron shower in the neutrino target, which is a 60 cm long aluminium 
cylinder. Although the final calculations were made using the parameterization, it 
is interesting to compare the GEANT yields with the ones measured in [2]. These 
measurements are rather detailed, they include secondary particles from the pri- 
mary beam of 67 - GeV protons in the momentum range from 5 - 7 GeV to 60 GeV 
and the angle range from to 50 mrad, subdivided into energy bins of 1 GeV and 
angle bins of 2 mrad. 

The comparison for GEANT 321 with GHEISHA is shown in Fig. 1. 

The comparison for GEANT 321 with FLUKA is shown in Fig. 2. 

From these fugures one can see that GEANT does not simulate correctly the 
yields of hadrons, the difference reaching a factor of 2 for the case of GHEISHA. 
For this reason in the mode of work with a full simulation the correction by weights 

(previous section, 4a) should be used. 

Special attention should be paid to the spectrum of K^, for which there are no 
reliable measurements at 70 GeV. Often it was assumed that there is the following 
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Fig. 1. The spectra of secondary particles from the thick aluminium target exposed to the 
beam of 67 - GeV protons, integral over the angle range 0-30 mrad. GEANT 321 with 
GHEISHA (CERN library 95a) - points with errors, line - measurements [2] 

relation between the charged and neutral K mesons: 

Kl = Kl=^ (1, 



It is clear from the following considerations that the relation should be different. 
We are usually interested in spectra of rather fast particles, which means that they 

should be "leading" particles, i. e. should contain a quark from the incident proton. 
Simple counting of quarks shows that the spectrum of should be closer to the 
spectrum of than to the spectrum of K^. 

In order to find the best description the following measurements were analysed: 
the measurement at FNAL [3] at 300 GeV, the charged kaon measurements at 
CERN at 400 GeV [4] and the charged kaon measurements [2] at IHEP (Protvino) 
at 67 GeV. All measurements were made on thin targets. The method [5] of compar- 
ison of spectra at different energies was used. It uses the variable xr = E*/E^^^, 
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Fig. 2. The spectra of secondary particles from the thick aluminium target exposed to the 
beam of 67 - GeV protons, integral over the angle range 0-30 mrad. GEANT 321 with 
GEANT - FLUKA (CERN library 95a) - points with errors, Une - measurements [2] 

the ratio of the energy in the center of mass system to the maximal energy at a given 
Pt. The best description was obtained with the following formula: 



where n is the ratio of the distributions u/d in a proton. This value is a function of 
a scaling variable x, which we substitute here by xji, which does not introduce a 
large error since the function does not change rapidly with x. It can be found for 
example in [6]. n{xR) changes from 1 to 5 as xr changes from to 1 (Fig. 3). The 
following approximation was used afterwards in calculations of spectra: 

n = -lA35x% + 5.3883a;^ + 0.22269a;R + 1.0189 (3) 



The uncertainty of formula (2) is estimated to be 15%, which is a quadratic sum 
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I The ratio u/d in a proton 



u/d 




Fig. 3. The ratio of the distributions u/d in a proton [6] 

of the uncertainty of the FNAL measurements (10%), the uncertainty of the scaling 
from 70 to 400 GeV (10%) and a smaller uncertainty of the CERN measurements. 

In Fig. 4 one can see that the formula (2) is in a much better agreement with 
the copper data than + K~)/2. The agreement of the formula (2) on berillium 
is also quite reasonable. Apparently the measurements [2] on berillium slightly 
underestimate the yields. 



4 The beam dump neutrino spectra. 

Practically the only possible way of spectra calculation for beam dump exper- 
iments is the mode with a full simulation of hadron showers in the target. It is 
difficult here to use parameterizations. However, at least the GEANT (GHEISHA) 
yields of particles from the first interaction of protons were compared with the re- 
sults of measurements on the thin copper target [2]. The difference was similar to 
that seen in Fig. 1. A corresponding correction was used in the calculations of spec- 
tra (see section 2, 4a). 

The systematic uncertainties of the calculated spectra were determined by the 
systematic uncertainties of the measurements [2]. Possible uncertaintied in the 
simulation of the secondary interactions by GEANT were suppressed since corre- 
sponding secondary particles contribute 20 - 30% to the neutrino spectra interesting 
for the beam dump experiment. These contributions were so small due to the fact 
that the experiment [7] was devoted to the search of hard neutrinos from the charm 
production which required rather high energy cuts. 

It is possible to reduce these uncertainties for the electron neutrino and antineu- 
trino spectra. For this purpose the reconstructed muon neutrino events were used. 
The high energy region of the muon neutrino and antineutrino spectra is dominated 
by decays of charged kaons (the contribution from charm estimated from the beam- 
dump extrapolation method was small). For this reason, comparing the number of 
reconstructed muon neutrino events with the calculated one it is possible to check 
the charged kaon spectra and to reduce systematic uncertainties. 
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Fig. 4. a) The line - spectrum measured at FNAL; circles - K'^ spectrum calculated 
according to the formula (2) from the charged kaon spectra measured at IHEP on a cop- 
per target; squares - spectrum calculated according to the traditional formula (1). b) 
The line - spectrum measured at FNAL; circles - spectrum calculated according 
to the formula (2) from the charged kaon spectra measured at IHEP on a berillium tar- 
get; squares - K° spectrum calculated according to the formula (2) from the charged kaon 
spectra measured at CERN on a berillium target. All spectra are at a fixed pt relevant for 
neutrino spectra. Uncertainties of neutral kaon spectra are shown, uncertaintied of charged 
kaon spectra are not shown. 

The final accuracy of the spectra calculations were 10% for the electron neutri- 
nos and 16% for the electron antineutrinos. 

Besides, for the beam-dump experiment such important value as a ratio of ef- 
fective densities of the two targets can be calculated with a very good accuracy. 
Since the dimensions of the half-density target were finite and there were construc- 
tive gaps, it is not exactly 2 as it was designed. Furthermore, this ratio for muon 
neutrinos is different from that for electron neutrinos since for the latter there is 
a contribution from K'^. For such calculations mainly a correct simulation of the 
shower shape is needed, which is reproduced by GEANT more reliably (this was 
checked, for example, during the calibration of the neutrino detector). The calcu- 
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lated ratio of the effective densities is 2.0 for muon neutrinos and 1.9 for electron 
neutrinos (accuracy 1%, mainly statistical). 

For additional control of the beam and spectra during the beam-dump experi- 
ment the muon fluxes in 9 gaps of the muon filter were measured [8]. In Fig. 5 the 
calculated fluxes are compared with the measured ones (the accuracy of measure- 
ments 3%, to which the uncertainty of the prompt muon pair production subtraction 
7% should be added). Unfortunately, such control does not allow to improve the ac- 
curacy of the electron neutrino spectra calculations since the main source of muons 
(more than 80%) is the decay of pions. 




Fig. 5. Beam-dump muon fluxes x 1000, integrals inside the circle with a radius 70 cm, per 
one incident proton. Line is the result of the calculation. 



5 The neutrino beam with the focusing system turned off. 

Two rather long runs of the IHEP-JINR neutrino detector were carried out with 
the neutrino channel focusing system turned off. Such mode of work has the advan- 
tage of a higher fraction of electron neutrinos, sufficient to perform some studies. 
These runs were also used for the measurement of the total cross section of neutrino 
and antineutrino interactions (simultaneously). The corresponding total neutrino 
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fluxes are significantly lower than with the focusing system turned on, but the loss 
of statistics for the neutrino detector was not dramatic because due to the lack of a 
resonant beam extraction (only a fast extraction is available in IHEP Protvino) the 
event aquisition rate is limited by the maximal allowed occupancy of the neutrino 
detector. In these measurements it is also important that neutrino and antineutrino 
spectra are less steep, which decreases the effects of smearing. 

Since in this neutrino beam (as well as in all the beams described below) the 
standard neutrino target was used, the calculations were accomplished in the mode 
of work of the program that uses the parameterization of secondary particle spectra 
from the target. But the data of [2] could not be used directly since the measure- 
ments [2] were carried out with a 6 cm diameter aluminuim target, whereas in the 
neutrino runs a 1 cm diameter target was used. The neutrino spectra in this case are 
better (harder), but the parameterization should be corrected. The correction was 
obtained in the following way: the program was run in the mode of full shower 
simulation and the correction was obtained by dividing the corresponding yields of 
hadrons from the targets. 

In the mode of work with the parameterization the GEANT correction proce- 
dure mentioned in section 2, 4a, is still needed since part of protons traverse the 
target without interaction and interact in the lenses, which is simulated by GEANT. 
What is significantly better in this mode of work is that statistical errors at high 
energy are much lower since the unweighted spectrum of hadrons from a target to 
be traced by GEANT is flat (is not a function of momentum). 

It turned out that without focusing the secondary effects, such as interactions of 
hadrons in the channel elements are much more important than with normal focus- 
ing. This happens because a bigger fraction of particles hit the walls of the decay 
tunnel. In addition, when the focusing system is turned on, the secondary particles 
from the interactions in the neutrino channel elements are not focused or are not in 
optimal conditions for a focusing. For this reason they contribute less to the neu- 
trino spectra. 

Only after calculations with a program based on GEANT the agreement be- 
tween the calculated and measured [8] muon fluxes was obtained. This is shown in 
Fig. 6. 

The systematic errors of the muon neutrino spectra in the first approximation 
are those of [2], but they are additionally controlled, mainly below 10 GeV (pionic 
part), by the muon flux measurements (the accuracy of the measurements about 

3%). 

In the electron neutrino spectrum for this neutrino channel a contribution from 
muon decays is significant. For the calculation of this contribution it is important 
to take into account the polarization of muons. These decays increase the electron 
neutrino spectrum by up to 20% at 3 GeV (thay give softer neutrinos so this differ- 
ence is the smaller the higher is the energy). This effect is even more important for 
the neutrino beam with focusing. 
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Fig. 6. Muon fluxes x 1000 for the neutrino beam with focusing system off. Solid line - 
calculations, circles - measurements, dashed line - calculations without secondary effects. 

6 The neutrino beam with the focusing system on. 



Here some difficulty was to introduce into the program the parabolic lenses used 
in IHEP because there is no paraboloid shape in GEANT. If define these lenses as 
cylinders, the magnetic field in them is not only non-uniform but non-contiguous 
and GEANT does not trace particles correctly. For this reason a routine for the step 
limitation in the lenses was introduced. The step is limited so that it is either en- 
tirely in a region with field or entirely in a region without field or is very small. 

The electron neutrino spectrum with a contribution from decays of muons is 
shown in Fig. 7. 

The systematic uncertainty of the neutrino flux calculation consists of the same 
uncertainties as for the beam with the focusing system off and additional uncer- 
taintied of the focusing system such as a possible inaccurate target and lenses po- 
sitioning and an uncertainty in the simulation of starting points of particles inside 
the target. No studies of the focusing system uncertainties were made with this 
program, but the previous studies [9] show that, taking into account muon flux 
measurements, they can be kept as low as 7% 
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Fig. 7. Vf, spectrum in the experiment with focusing. Dashed Une - the same without taking 
into account decays of muons. 



7 The neutrino beam with a short decay base. 



For this experiment [10], with rather complicated and asymmetrical geometry 
of the decay region, it was very convenient to use GEANT tools for the description 
of the geometry. 

The decay cavity can be described as 4 adjacent (with steps in the vertical 
projection) trapezoids, asymmetric in the vertical projection. It is narrow near the 
target to provide a better radiation shielding and wider near the end closest to the 
detector. It was made asymmetric because of the danger of "underground" muons 
that can go around the iron muon filter, in the earth under it due to a multiple 
scattering. 

In this experiment a contribution from muon decays is also rather significant, 
including decays of muons in the muon filter (the path length of muons in it is 
comparable or even bigger than the length of the decay region). 
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8 Conclusion 

A program for the neutrino spectra calculations based on the GEANT library is 
created. The range of possible applications of the program is very wide. 
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